A risk model based on 10 ferroptosis regulators and markers established by LASSO-regularized linear Cox regression has a good prognostic value for ovarian cancer patients

Ovarian cancer is the deadliest gynecologic cancer due to its high rate of recurrence and limited early diagnosis. For certain patients, particularly those with recurring disorders, standard treatment alone is insufficient in the majority of cases. Ferroptosis, an iron- and ROS (reactive oxygen species)-reliant cell death, plays a vital role in the occurrence of ovarian cancer. Herein, subjects from TCGA-OV were calculated for immune scores using the ESTIMATE algorithm and assigned into high- (N = 185) or low-immune (N = 193) score groups; 259 ferroptosis regulators and markers were analyzed for expression, and 64 were significantly differentially expressed between two groups. These 64 differentially expressed genes were applied for LASSO-regularized linear Cox regression for establishing ferroptosis regulators and a markers-based risk model, and a 10-gene signature was established. The ROC curve indicated that the risk score-based curve showed satisfactory predictive efficiency. Univariate and multivariate Cox risk regression analyses showed that age and risk score were risk factors for ovarian cancer patients’ overall survival; patients in the high-risk score group obtained lower immune scores. The Nomogram analysis indicated that the model has a good prognostic performance. GO functional enrichment annotation confirmed again the involvement of these 10 genes in ferroptosis and immune activities. TIMER online analysis showed that risk factors and immune cells were significantly correlated. In conclusion, the risk model based on 10 ferroptosis regulators and markers has a good prognostic value for ovarian cancer patients.


Introduction
Ovarian cancer is the deadliest gynecologic cancer due to its high rate of recurrence and limited early diagnosis [1], putting a huge cost upon patients and society.There has been a considerable decrease in the incidence and death rates of ovarian cancer over the last several decades because of advancements in therapy; however, the death rate from ovarian cancer remains high, and less than one-half of patients survive for more than five years [2].Modern treatments for ovarian cancer patients vary with individual specialty, depending on histological type and cancer stage.In general, the standard of care for patients with ovarian cancer consists of primary debulking surgery and platinum-based combination chemotherapy [3].Nevertheless, for certain patients, particularly those with recurring disorders, standard treatment alone is insufficient in the majority of cases.
Ferroptosis, a distinct form of programmed cell death, was first described in 2012 by Dixon [4].Ferroptosis, in contrast to autophagy and apoptosis, is defined as an iron-and ROS (reactive oxygen species)-dependent cell death characterized primarily by cytological alterations, such as diminished or disappeared mitochondrial cristae, rupture of outer mitochondrial membrane [5][6][7][8][9].These cellular abnormalities were caused by a loss of plasma membrane selective permeability as a result of severe membrane lipid peroxidation and the initiation of oxidative stress [10].Ferroptosis process regulation involves various genes, which could be briefly grouped into drivers that drive ferroptosis, suppressors that suppress ferroptosis, and markers that indicate ferroptosis occurrence [11].Accumulating evidence has indicated that ferroptosis is closely correlated with various diseases, including neurodegenerative disorders [12], ischemia/ reperfusion damage [13,14], acute renal injury [15], and malignancies, etc.Thus, ferroptosis regulators and markers might be potential targets and markers for diseases, including cancers.
A large body of data suggests the inhibitory effect of ferroptotic cell death on tumor development.Despite the fact that oxidative phosphorylation is most effective in the production of ATP, numerous tumor cells undergoing metabolic reprogramming primarily produce ATP from cytosolic aerobic glycolysis combined with lactate fermentation.This metabolic reprogramming in malignancies was notably found in the 1920s by Warburg and Cori and has been proposed as a tumor cell way to prevent toxic ROS levels [16,17].Nevertheless, maintaining this Warburg effect necessitates elevated glucose uptake and increased metabolic activity which make cancer cells significantly dependent on the anti-oxidant mechanism and possibly even more vulnerable to oxidative stress [18,19].As a result, highly proliferative tumor cells have been found to require handling of increased ROS levels in order to effectively develop tumors [20][21][22].Thus, cancer cells have a much higher demand for iron than noncancer cells and such reliance upon iron will make tumor cells more sensitive to ferroptosis.Cancer cells also have ROS tolerance and dependence on iron in metabolism, making them more susceptible to ferroptosis.Given these previous findings, ferroptosis regulators (drivers and suppressors) and markers might be promising and therapeutic targets for cancer.Regarding ovarian cancer, a potential link between ferroptosis and ovarian cancer is demonstrated based on previous studies [23,24], and a focus on ferroptosis may provide a sophisticated therapeutic strategy for treating ovarian cancer [25].Immunotherapy, including immune checkpoint inhibitors, has been used to treat ovarian cancer [26].It is clinically important to investigate the factors affecting the prognosis of immunotherapy.Accumulating evidence suggests that ferroptosis plays a crucial role in immune evasion [26].In turn, in immunotherapy, CD8 + T cells can trigger ferroptosis in cancer cells [27].Targeting ferroptosis in combination with immunotherapy might become a prospective strategy for cancer therapy.
Although several studies have proposed a variety of prognostic models based on differentially expressed genes between ovarian cancer and para-cancerous samples, the sensitivity and specificity of these prediction models remain unsatisfactory.Recently, the Least Absolute Shrinkage and Selection Operator (LASSO) has been proposed as a regression algorithm for highdimensional data [28], which was applied to select the most prominent predictive characteristics in the training dataset.Herein, cases of ovarian cancer patients with clinical information and expression profile information were obtained from The Cancer Genome Atlas Ovarian Cancer (TCGA-OV).Cases were calculated for immune scores using the Estimation of Stromal and Immune cells in Malignant Tumors using Expression data (ESTIMATE) algorithm [29], ferroptosis regulators and markers were obtained from the FerrDb database (http:// www.zhoun an.org/ ferrdb/), and differentially expressed ferroptosis regulators and markers between high-and low-immune score groups were analyzed.Then, ferroptosis regulators and markers associated with ovarian cancer patients' OS were identified using the univariate and multivariate Cox regression analyses.Next, The TCGA-OV cases were randomly divided into a training set and a validation set.The LASSO-regularized linear Cox regression was employed to construct a risk model consisting of ferroptosis regulators and markers.The risk model's prognostic value was then validated using a time-dependent receiver-operating characteristic (ROC) curve analysis, multivariate Cox's proportional hazard regression model analysis, and nomograms.The LASSO-regularized linear Cox regression was used to establish a ferroptosis regulator-and markers-based risk model for ovarian cancer prognosis.

Source data of OV from TCGA and GEO database
All data were collected from the Cancer Genome Atlas (TCGA) database (https:// tcga-data.nci.nih.gov/ tcga/), including gene expression profiles and clinical data of patients enrolled.The inclusion criteria of the present study included: (1) Patients with complete gene expression profiles that could be applied to assess immune scores using the ESTIMATE algorithm; (2) patients with complete clinical data; and (3) patients with complete prognosis data.Then, a total of 378 cases of ovarian cancer patients with clinical information and expression profile information were obtained from TCGA-OV database.
The microarray data have been deposited in NCBI's Gene Expression Omnibus and are accessible through Gene Expression Omnibus (GEO; https:// www.ncbi.nlm.nih.gov/ geo/) Series accession number.Gene expression data have been archived as microarray datasets in GEO repository at the NCBI archives and are accessible through GEO Series accession number GSE63885.The GSE63885 dataset includes the gene expression profiles from 101 ovarian carcinoma specimens: 73 serous, 12 endometrioid, 9 clear cells, and 7 undifferentiated.

Data processing
A total of 378 cases of ovarian carcinoma patients with clinical information and expression profile information were obtained from The Cancer Genome Atlas Ovarian Cancer (TCGA-OV).Each patient's immune score was evaluated by the ESTIMATE algorithm [29].According to this score and using the mean score as a cut-off, ovarian cancer patients were assigned into high immune score (N = 185) and low immune score (N = 193) groups for Kaplan-Meier survival analysis.Subsequently, based on machine learning, these 378 patients with ovarian cancer were randomly divided into a training set (N = 189) and a validation set (N = 189).The clinical information is shown in Table 1.

Identification of prognostic signature
The LASSO-regularized linear Cox regression was implemented using the sklearn library in Python for the Lasso regression model.For univariate and multivariate Cox regression analyses, we utilized the survival package in the R language.Based on 189 samples of ovarian cancer patients in the training set, the prognostic value of differentially expressed ferroptosis regulators and markers was evaluated by the multivariate Cox regression analyses identifying ferroptosis regulators and markers associated with ovarian cancer patients' overall survival.Then, the coefficients of these ferroptosis regulators and markers linked to ovarian cancer patients' overall survival were analyzed using the LASSOregularized linear Cox regression.The following formula was employed to calculate each patient's risk score for prognostic signature: risk score = expression of gene 1 × β 1 gene 1 + expression of gene 2 × β 2 gene 2 + expression of gene 3 × β 3 gene 3 [30,31].For the present study, the formula is: Risk_score = Exp(LAMP2)*-

Genomic-clinicopathologic nomogram
The nomogram was constructed based on the LASSOregularized Cox regression model using the "survival" and "rms" package in the R language.The Consistency index (C index) was calculated using the survConcordance function in "Survival" package in R language, which is the capability of the model to distinguish between patients who survived and those who did not.The performance of the nomogram was evaluated using a bootstrap resampling approach with 1000 iterations to estimate the bias-corrected C-index.The bootstrap samples were generated by randomly sampling with replacement from the original dataset, using the same sample size as the original.This process was repeated 1000 times to ensure the stability and reliability of the nomogram performance estimates.Then, the calibration curves for the probability of 1-, 3-. and 5-year overall survival (OS) showed satisfactory agreement between predicted survival and actual observed survival.

Sample collection
Ovarian cancer samples were collected from 20 ovarian carcinoma patients in Affiliated Tumor Hospital, Xinjiang Medical University.The inclusion criteria for ovarian cancer patients comprise: (1)  prior to surgery; (5) no history of other malignancies or ovary-associated disorders.The normal ovarian tissue specimens were obtained from 10 cases who received myomectomy.All tissue samples were fixed in formalin until use.Written informed consent was obtained from all patients enrolled.Human tissue experiments were approved by the Ethics Committee of Affiliated Tumor Hospital, Xinjiang Medical University (approval ID: K-2021057).The clinicopathologic characteristics of 10 low-grade (pathological grade G1-G2) and 10 high-grade (pathological grade G3) ovarian carcinoma patients were included in Table 7.

Statistical analyses
All statistical analysis uses Python (version 3.7.5;https:// www.python.org/) and R language (version 4.0.2;https:// www.r-proje ct.org/).P < 0.05 is considered statistically significant.Delete the missing clinical data from the list; delete the entire sample from the analysis if the value of any parameter is missing.OS is defined as the time interval between the date of the first patient visit and the date  of death.The mean value comparison of continuous variables uses a two-sided t-test.Kaplan-Meier method was used for survival analysis of high-risk group and low-risk group using a two-sided log-rank test in Python.The log-rank test, a non-parametric test used to compare the survival distributions of two groups, was employed to assess the significance between high and low-risk groups in the Kaplan-Meier survival analysis.This test is particularly suited for censored data and is widely used in survival analysis to test the null hypothesis that there is no difference in survival between the groups being compared.Patients with risk scores above the median were classified into the high-risk group, whereas those with scores below the median were classified into the low-risk group.We conducted two-sided t-tests for the hypothesis testing of 259 comparisons.This choice was made to account for the possibility of both positive and negative effects.FDR correction was performed with a cut-off value of 0.05 [32].

Ferroptosis regulator and marker-based risk model established using the LASSO-regularized linear Cox regression
First, a total of 378 ovarian cancer subjects with clinical information and expression profile information from TCGA-OV were calculated for immune score with ESTI-MATE algorithm and assigned into high-or low-immune score group (Table 1) using the mean immune score as the cut-off; a Kaplan-Meier estimate was employed to analyze the correlation between immune score and ovarian carcinoma patients' OS.As shown by Fig. 1A, patients from the high-immune score group exhibited an obviously better prognosis.
Then, the list of the ferroptosis-related regulators and markers was retrieved from the FerrDb database (http:// www.zhoun an.org/ ferrdb/).A total of 259 ferroptosisassociated regulators and markers were analyzed for differential expression in the high-or low-immune score group and 64 were found with significant differences  2).Based on 189 ovarian cancer samples in the training set, the prognostic value of 64 ferroptosis-related regulators and markers was evaluated by multivariate risk regression analysis, and 15 of them showed to be remarkably linked to ovarian carcinoma patients' OS (P < 0.05) (Table 3).
Subsequently, 378 ovarian cancer patients from TCGA-OV were randomly separated into a training set (N = 189) and a validation set (N = 189) (Table 1).For predicting the clinical outcome of these ferroptosis-related regulators and markers, a LASSO-regularized linear Cox regression analysis was carried out to establish the risk score model based on 189 samples from the training set.According to the minimum criteria, a 10-gene risk signature was established (Fig. 1C).The 10 ferroptosis-related regulators and markers are: LAMP2, NOS2, ALOX5, CD44, CHMP5, FH, GOT1, DUOX2, SLC7A11, and DDIT3.

The prognostic value of the risk score model
According to the regression coefficient, the risk score for each subject in the training set, validation set, and independent dataset GSE63885 was calculated with the aforementioned formula.Subjects in each dataset were then assigned into high-and low-risk score groups based on the median score in each group.In the TCGA-OV training set (Fig. 2A), TCGA-OV validation set (Fig. 2B) and TCGA-OV entire set (Fig. 2C), and independent dataset GSE63885 (Fig. 2D), the subjects with lower risk scores obtained better overall survival.For further confirming the prognostic value of the risk

Univariate and multivariate Cox analyses of clinicopathological variables
Next, the clinical characteristics in TCGA-OV patients were analyzed using a univariate and multivariate Cox's proportional hazard regression model (Table 4).Risk scores were considerably correlated with OS in OA patients.Also, Fig. 3A showed that according to TCGA-OV data, age (P < 0.001, HR = 1.023; 95% CI = 1.010-1.04)and the risk score (P < 0.001, HR = 1.029; 95% CI = 1.017-1.04)could predict patients' OS.Moreover, patients from TCGA-OV were assigned into high-risk score or low-risk score group and the immune scores were shown.Figure 3B showed that patients in the high-risk score group obtained lower immune scores.
For validating the prognostic value of the risk score model, a Nomogram analysis was performed based on TCGA-OV. Figure 4A-D showed that based on the cases with prognostic information from TCGA-OV, a prognostic nomogram predicting the 1-, 3-, and 5-year survival probability was established, respectively.The nomogram included age, stage, and the risk score.The C-index of the risk score was 0.65, suggesting a favorable prognosis of the model.

Gene Ontology (GO) analysis on the 10 ferroptosis regulators and markers
For an in-depth understanding of the 10 genes forming the risk model, GO analysis was conducted.Metascape online tool (https:// metas cape.org) [34] confirmed again the involvement of these 10 genes in ferroptosis such as cell metabolism, oxidation-reduction reactions, and oxidative stress, as well as immune activities such as leukocyte activation, granulocyte activation, and neutrophil-mediated immunity (Fig. 5; Table 5).
Since the GO analysis confirmed again that these 10 genes were enriched in ferroptosis biological activities and immune activities, next, the association of 10 risk factors with immune cells was analyzed in ovarian cancer.Through TIMER online analysis, the association of immune cells with ovarian cancer patients' cumulative survival was analyzed using a Kaplan-Meier estimate; Fig. 6A showed that patients with higher dendritic cell or CD4 + T cell obtained better cumulative survival.Similarly, patients with higher DDIT3 or CHMP5 obtained better cumulative survival (Fig. 6B).Pearson correlation coefficient analysis was performed to analyze the correlation between immune cells and the risk factors; Fig. 6C; Table 6 showed that risk factors and immune cells were significantly correlated.

Expression of immune microenvironment and ferroptosis markers in ovarian cancer and its relationship with prognosis of ovarian cancer patients
To validate the bioinformatics analysis results, IHC staining was applied to investigate the expression of immune microenvironment markers (dendritic cell CD CD1α, CD4 T cell, CD4) and ferroptosis markers (CHMP5 and DDIT3) in normal ovarian tissues and low-or high-grade ovarian cancer tissues (Fig. 7).As showed in Fig. 7A-B, compared to normal ovarian tissues, CD1α and CD4 protein was markedly low-expressed in low-or high-grade ovarian cancer tissues; and CD1α was even lower in high-grade ovarian cancer tissues.For ferroptosis markers (CHMP5 and DDIT3), CHMP5 (Fig. 7C) and DDIT3 (Fig. 7D) protein expressions were notably down-regulated in high-grade ovarian cancer tissues.The clinicopathologic characteristics of clinical samples were shown in Table 7.

Discussion
Herein, subjects from TCGA-OV were calculated for immune scores using the ESTIMATE algorithm and assigned into high-(N = 185) or low-immune (N = 193) score group; 259 ferroptosis regulators and markers were analyzed for expression and 64 were significantly  10 genes in ferroptosis and immune activities.Pearson's correlation analysis showed that risk factors and immune cells were significantly correlated.
A growing body of data suggests that innate and adaptive immune systems play a critical role in the occurrence and progression of malignancies [35][36][37].Furthermore, cancer immunotherapy has progressed so rapidly that immune-based ovarian cancer prognostic signatures might provide potential value for identifying new molecular targets [38].By using the ESTIMATE algorithm, cases from TCGA-OV were calculated for immune scores and assigned into high-and low-immune score groups, and a Kaplan-Meier estimate then indicated that patients with higher immune scores obtained better overall survival, suggesting the prognostic value of the immune score.Given the crucial role of ferroptosis in ovarian cancer [39], 259 ferroptosis regulators and markers were obtained from the FerrDb database (http:// www.zhoun an.org/ ferrdb/) [11] and differentially expressed ferroptosis regulators and markers between the immune score high group and low group were evaluated.Between the two groups, a total of 64 ferroptosis regulators and markers were differentially expressed, suggesting those ferroptosis-related genes are associated with the level of infiltrating of stromal and immune cells in ovarian cancer and might consist of a signature for ovarian cancer prognosis.Consistent with the hypothesis, interferon-γ has been reported to induce ferroptosis within tumor cells [40].Moreover, the early ferroptotic cancer cells might significantly promote immune responses [41].
Exploiting ferroptosis inducers provides a potential therapeutic method for the treatment of ovarian cancer.Multiple conventional drugs might trigger ferroptosis in tumor cells, including Sulfasalazine, Artesunate, Temozolomide, and Cisplatin [42].Similarly, ferroptosis suppressors represent promising therapeutic targets for treating ovarian cancer.A combination of Ferroptosis inducers combined with chemotherapeutic agents gains a remarkable synergistic effect on their anti-tumor activity [43].
In the present study, 15 of the 64 ferroptosis regulators and markers were considerably correlated with ovarian carcinoma patients' OS.Furthermore, the LASSO-regularized linear Cox regression established a 10-gene risk model predicting ovarian cancer prognosis, consisting of LAMP2, NOS2, ALOX5, CD44, CHMP5, FH, GOT1, DUOX2, SLC7A11, and DDIT3, all of which could play a role in ovarian cancer development.For instance, cardamonin inhibited the mTOR lysosomal colocalization LAMP2, a well-known protein found in the membrane of lysosomes, decreasing Raptor siRNA SKOV3 cell proliferation [44].NOS2 has initially been reported to play a significant anti-tumor role in the immune response; however, growing evidence has demonstrated that the expression level of NOS2 in tumor cells is usually associated with impaired prognosis [45].Inherited variation in ALOX5 seems to affect ovarian cancer risk [46].Exosomal CD44 has been reported to enhance the capacity of ovarian cancer cells to invade via CD44 transfer to the peritoneal mesothelium [47].GOT1 modulates cellular metabolism by coordinating the use of carbohydrates and amino acids to satisfy dietary needs for long-term proliferation; in ovarian cancer, adapalene suppresses the growth of ovarian cancer cells by binding to GOT1 [48].PARP pharmacologic inhibition represents a major factor in treating ovarian cancer with mutations in the BRCA; inhibition of PARP enhances ferroptosis through suppressing SLC7A11 and synergizes with ferroptosis inducers within BRCA-proficient ovarian cancer [49].By grouping TCGA-OV training set and validation set and GSE63885 samples into 2 sub-groups with the median value of risk score as a cut-off, the Kaplan-Meier survival estimate indicated that the risk score was strongly associated with the overall survival of patients, the ROC curve demonstrated that the risk score-based curve showed satisfactory prediction efficiency, and the multivariate Cox's proportional hazard regression analysis identified the risk score as an independent risk factor.More importantly, the intuitive and effective nomogram integrating age, stage, and risk model was established, and it could be used to predict the outcomes of patients quickly.Through TIMER online analysis and Kaplan-Meier analysis, patients with higher levels of dendritic cell, CD4 + T cell, DDIT3, and CHMP5 were associated with better cumulative survival.The clinical samples also confirmed that the levels of immune microenvironment markers (CD1α and CD4) and ferroptosis markers (CHMP5 and DDIT3) were lower in high-grade ovarian cancer tissues.These mRNAs might be used as therapeutic targets in treating ovarian cancer.
Regarding the limits associated with this research, firstly, the biological effect of the 10 identified genes should be verified by in vitro and in vivo experiments; considering the 10 genes were ferroptosis regulators and markers Taken together, it is believed that the risk model based on 10 ferroptosis regulators and markers has a good prognostic value for ovarian cancer patients.It is worth noting that the risk score can also significantly distinguish ovarian cancer from normal samples, which may have a certain auxiliary value for early clinical screening of ovarian cancer.• thorough peer review by experienced researchers in your field • rapid publication on acceptance • support for research data, including large and complex data types • gold Open Access which fosters wider collaboration and increased citations maximum visibility for your research: over 100M website views per year

•
At BMC, research is always in progress.

Learn more biomedcentral.com/submissions
Ready to submit your research Ready to submit your research ?Choose BMC and benefit from: ?Choose BMC and benefit from:

Fig. 1
Fig. 1 Ferroptosis regulator and marker-based risk model established using the Least Absolute Shrinkage and Selection Operator (LASSO)-regularized linear Cox regression.A A total of 378 patients were calculated for immune score by ESTIMATE algorithm and assigned into high-or low-immune score group; the association of immune score with the overall survival in patients with ovarian cancer was analyzed using a Kaplan-Meier estimate.B Ferroptosis regulators and markers based on FerrDb database (http:// www.zhoun an.org/ ferrdb/) were analyzed for differential expression in high-or low-immune score groups.C The coefficients of differentially expressed ferroptosis regulators and markers in TCGA-OV training set (n = 189) were calculated by multivariate Cox regression using LASSO-regularized linear Cox regression

Fig. 2
Fig.2The prognostic value of the risk score model.A Risk scores of cases in TCGA-OV (training set, N = 189) were calculated using the formula described in the M&M section.Cases were divided into two groups using the median value of the risk score as a cut-off; the correlation of the subjects' overall survival with risk score was analyzed.B Risk scores of cases in TCGA-OV (validation set, N = 189) were calculated using formula described in the M&M section.Cases were divided into two groups using the median value of the risk score as a cut-off; the correlation of the subjects' overall survival with risk score was analyzed.C Risk scores of cases in TCGA-OV (all subjects, N = 378) were calculated using formula described in the M&M section.Cases were divided into two groups using the median value of the risk score as a cut-off; the correlation of the subjects' overall survival with risk score was analyzed.D Risk scores of cases in GSE63885 were calculated using formula described in the M&M section.Cases were divided into two groups using the median value of the risk score as a cut-off; the correlation of the subjects' overall survival with risk score was analyzed.E Receiver operating characteristic (ROC) curves showed the predictive efficiency of the risk score model in 1-, 3-, 5-, 8-, or 10-year overall survival based on cases from TCGA-OV.F ROC curves showed the predictive efficiency of the risk score model in 1-, 2-, 3-, 4-, or 5-year overall survival based on cases from GSE63885

Fig. 3
Fig. 3 Cox multivariate analyses of clinicopathological variables.A The multivariate Cox regression analyzing the association between ovarian cancer patients' overall survival and age, stage, and risk score based on TCGA-OV data.B Patients from TCGA-OV were assigned into high-risk or low-risk group and the immune scores were shown

Fig. 4
Fig. 4 Nomogram analysis based on TCGA-OV.A Nomogram composed of age, stage, and risk score for the prediction of 1-, 3-, and 5-years survival probability.B-D Calibration plot for the evaluation of the nomogram in predicting 1-year, 3-years, and 5-years overall survival

Fig. 6
Fig. 6 Association of 10 risk factors with immune cells in ovarian cancer.A The association of immune cells with ovarian cancer patients' cumulative survival was analyzed using a Kaplan-Meier estimate.B The association of each of the 10 risk factors with ovarian cancer patients' cumulative survival was analyzed using a Kaplan-Meier estimate.C The correlation between immune cells and the risk factors was analyzed using Pearson correlation coefficient.The left panel: heatmap showing the significance of the correlations (P-value); the right panel: heatmap showing the strength and direction of the correlations (R-value)

Table 1
Clinical information of patients in the high-, low-immune score groups, training set, and validating set diagnosed as serous ovarian cancer by postoperative pathology; (2) > 18 years old; (3) received ovarian cancer resection; (4) no chemotherapy, radiation therapy, immunotherapy, and other therapies

Table 2
Differentially-expressed ferroptosis regulators in highand low-immune score groups

Table 3
Ferroptosis regulators associated with ovarian cancer patients' overall survival

Table 4
Cox risk regression analysis on the correlation between the overall survival of ovarian cancer patients and clinical features

Table 5
Gene ontology functional enrichment annotation differentially expressed between two groups.These 64 differentially expressed genes were applied for LASSOregularized linear Cox regression for establishing ferroptosis regulators-and markers-based risk model, and a 10-gene signature was established.The ROC curve indicated that the risk score-based curve showed satisfactory predictive efficiency.Based on univariate and multivariate Cox risk regression analyses, age, and risk score were risk factors for ovarian cancer patients' OS; patients in the high-risk score group obtained lower immune scores.The Nomogram analysis indicated that the prognostic outcomes of the model were consistent with the actual outcomes.GO functional enrichment annotation confirmed again the involvement of these

Table 6
Association of 10 risk factors with immune cells in ovarian cancer

Table 7
The clinicopathologic characteristics of normal control and ovarian cancer patients N Normal, P Ovarian cancer patient, FIGO International Federation of Gynecology and Obstetrics